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ABSTRACT 

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the 
entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with 
spatial resolution of ^140 pc and find that a variable dust-emissivity index (/3) is required to fit the 
data. We find no significant long-wavelength excess above this model suggesting there is no cold dust 
component. We show that the gas-to-dust ratio varies radially, increasing from ^20 in the center to 
^70 in the star-forming ring at lOkpc, consistent with the metallicity gradient. In the lOkpc ring 
the average j3 is ~1.9, in good agreement with values determined for the Milky Way (MW). However, 
in contrast to the MW, we find significant radial variations in f3, which increases from 1.9 at lOkpc 
to ^2.5 at a radius of 3.1 kpc and then decreases to 1.7 in the center. The dust temperature is fairly 
constant in the lOkpc ring (ranging from 17-20 K), but increases strongly in the bulge to ^30 K. 
Within 3.1 kpc we find the dust temperature is highly correlated with the 3.6 /im flux, suggesting 
the general stellar population in the bulge is the dominant source of dust heating there. At larger 
radii, there is a weak correlation between the star formation rate and dust temperature. We find 
no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an 
estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value 
of (1.9 ±0.4) x KPcm-^Kkms- 1 ]- 1 . 
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Astronomy at long infrared wavelengths (20-1000 /im) 
is a relatively young field due to the need for space mis- 
sions to avoid the absorption of the atmosphere in this 
waveband. This waveband, however, is vital for astro- 
nomical studies as this is where dust in the interstel- 
lar medium (ISM) radiates. This is important for stud- 
ies of galaxy evolution as star formation regions tend 
to be dusty, and therefore the use of UV and opti- 
cal measurements to trace the star fo rmation ra t e can 
lead to it being underestimate d fsee iKennicuttl Il998| ; 
plain et all [1991 [Caliettil [20011 IPapovich fe"Bel1 [20021 
ICalzetti et al.l 120101) . Calibrating the relationship be- 
tween infrared emission and star formation rate has been 
difficult due to uncertainties from the contribution of 
the general stellar population to heating the dust, the 
fact that not all optical/UV emission is absorbed, and 
uncertainties in the properties of the dust. The dust 
emission is affected by the composition of the dust, and 
the proportion of non-equilibrium to equilibrium heat- 
ing. Studies with previous space missions IRAS, ISO, 
Spitz er and AKARI tried to address these questions 
(e.g-JWalterbos fc Greenawaltl | 1996l iBoselli et al.ll2004l: 
ICalzetti et al.ll2010l : iBuat et al.ll2011[) . However, as they 
were limited to wavelengths less than 160 /xm, they were 
not sensitive to the cold dust (< 15 K) and missed up to 
50% of the dust mass in galaxies. 

The continuum emission from the dust has been 

* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 
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proposed as a pot ential method of measuring the total 
mass of the ISM (IHikiebrandl fl983l: IGuelin et al.l H993I: 
lEales et al.ll2010l 120121 ); traditionally the amount of gas 
has been measured using Hi and CO surveys, but due to 
sensitivity and resolution issues this method is limited to 
low re dshift and small numbers of galaxies. Smi th et al.1 
(|2012f ) found that for early- type galaxies (E/SO) the 
ISM was detected for 50% of objects through its dust 
emission but only 22% through its CO emission. In 
addition, the conversion of the CO tracer to molecular 
gas mass is highly u ncertain and is believed to vary 
(see iWisonl 119951 iBoselli et all [2 002; 
" 120051: iNaravanan et al.N2012D . 



Ilsraell 



with metallicity 
iStrong et all 120(1 
This is a topical area as recent studies using the Planck 
all-sky survey and IRAS m aps have made an estimate o f 
the amount of 'dark gas' (jPlanc k Colla boration! l2011a| ) 
in the Milky Way. The 'dark gas' is thought to be 
molecular gas which is traced by dust, but not detected 
with the standard CO method. Previous works have 
also suggested the presence of 'dark gas' by using 7-ray 
emission from cosmic-ra y interactio ns with clouds of 
gas (jGrenier et al.l 120051 : [Abdo et al.l [20T0I) and bv the 
kinem atics of recycled dwarf galaxies |Bournaud et al.l 
12001 . 

Herschel is one of the European Space Agency's flag- 
ship observatories and observ es in the far-infrared (FIR) 
in the range 55-671 /jm fsee iPilbratt et al.|[2010h . Due 
to the large space mirror and cryogenic instruments, it 
has a high sensitivity and unprecedented angular reso- 
lution at these wavelengths. Herschel has the ability to 
target both large numbers of galaxies and map large ar- 
eas of sky. It has two photometric instruments: PACS 
(|Poglitsch et al.l 12010ft which can observe in 3 broad 
bands around 70, 100 and 160 /im (70 and 100 urn cannot 
be used simultaneously) and SPIRE (Gri ffin et al.ll20iol ) 
which observes simultaneously in filter bands centered 
at 250, 350 and 500 /im. SPIRE provides flux measure- 
ments on the longer wavelength side of the FIR dust peak 
(^160 /im), allowing us to obtain a full c ensus of dust in 
nearby galaxies (e.g. JDunne et al.1l2011| ). 

Andromeda (M31) and the Milky Way (MW) are the 
only two large spirals in the Local Group. Studies 
of Andromeda therefore provide the best comparison 
to observations of our own Galaxy with the advantage 
that we get a 'global' picture, whereas investigations 
of the Milky Way are limited by our location within 
the Galaxy. The total size of M31 and the scalelength 
of its dis k are both ap proximately twice those of the 
MW (see I Yin et al.l [20091 and references therein). How- 
ever, the star-formation rat e of the MW is ~3-6 M Q yr _1 
([Boissier fe Prantzosl 119991) c ompared to on ly ~0. 3-1.0 
in M31 (jBarmbv et al.l 120061: iWilliamsl 120031 G. P. Ford 
et al., i n preparation) , despite similar amounts of gas 
present (|Yin et al.l 120091) . For this reason, M31 is often 
labelled as 'quiescent'. The dust emission from M31 is 
dominated by a dusty star-forming ring at a radius of 
lOkpc, and was firs t observed in the infrared by IRAS 
([Habing et al.lU98l . 

Many projects to map the ISM in M31, have been un- 
dertaken using observations i n the mid-infrared (MIR) 

" 120061 ) . in the FIR with 
the Hi atomic line 



with Spitzer (B armby et al 



Spitzer dGordon et a 



2006) 



(|Thilker et al.l 120041 : iBraun et al.l ifM IChemin et all 



120091) and in the CO(J=1-0) line (INieten et al.1 12006T) . 
Studies of the gas kinematics and dust emission show 
M31 h as a complica t ed st ruc ture (e.g.. IChemin et al.l 
120091 : IGordon et all 120061 ): IBlock et al.l (|2006D at- 
tribute many of the features observed to density 
waves from a possible head-on collision with M32. 
iTabatabaei fc B crkhuiisenl (|2010l) investigated the rela- 
tion between ga s, dust and star fo rmation using FIR 
Spitzer data. iLerov et"aH ([201 ID used Spitzer and 
gas observations to investigate the conversion between 
CO(J=1-0) line flux and molecular hydrogen column- 
density in a sample of nearby galaxies including M31. 
In their analysis they found M31 has the lowest dust 
temperatures in their sample and therefore would bene- 
fit most by including Herschel data. 

The Herschel Exploitation of Local Galaxy Andromeda 
(HELGA) is a survey covering a ~ 5.5° x 2.5° area cen- 
tered on M31. We use PACS-SPIRE parallel mode, ob- 
serving at 100, 160, 250, 350 and 500 /im simultaneously. 
Further details of the observatio ns can be found in Sec- 
tion HH] and in lFritz et all ([201 ID . HELGA observations 
have been used in other work s to investigate s tructures 
in dust and Hi at large radii ([Fritz et al.l 1201 lh . the re- 
lationship between gas and star formation (G. P. Ford 
et al., in preparation), and the structure of M31 and the 
cloud-mass function (J. Kirk et al., in preparation). 

In this paper we use the Herschel data combined with 
the wealth of ancillary data to investigate the distri- 
bution, emission properties and the processes heating 
the dust in M31 on spatial scales of ~140kpc. There 
have still been relatively few attempts to map the dust 
within a ga l axy; r e cent studies wi t h Her s chel include 
iSmith et al.l (120101) iBoquien et ail (1201 ID. iFovle et al.1 
([2012D . IBendo et al.l (f2012t ). and lMentuch et al.l ([2012D . 
but they are often limited to small numbers of indepen- 
dent pixels. The advantage of M31 over these other stud- 
ies is the close proximity which allows us to investigate 
the dust at higher spatial scale and with many more in- 
dependent pixels. We also apply the Planck me thod for 
finding 'dark gas' ([Plan ck Collabora tion! l2011al ) to M31 
and use this method to determine a value of the X-factor 
— the relationship between the molecular hydrogen col- 
umn density and the observed CO tracer. In Section [2] 
we present the data used for this analysis and Section [3] 
describes our method for fitting the spectral energy dis- 
tribution of dust. Section [4] presents our maps of the 
dust properties, including the distribution of dust sur- 
face density, temperature and spectral index. This sec- 
tion also describes a comparison of the distributions of 
gas and dust and of the search for an excess emission 
at long wavelengths. In Section [5] we discuss the dust 
properties including the composition of the dust and the 
processes heating the dust. In this section we also search 
for 'dark gas' and determine the value of the X-factor. 
The conclusions are presented in Section [5] 

2. THE DATA 

2.1. Far-Infrared Observations 

Herschel observations of M31 were taken using both 
PACS and SPIRE in parallel mode covering an area of 
~ 5.5° x 2.5° centered on M31. To observe an area this 
large, the observations were split into two halves with a 
cross-scan on each half, which produced data at 100, 160, 
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250, 350 and 500 /im simultaneously (observation IDs: 
1342211294, 1342211309, 1342211319, 1342213207). Full 
details of t he observing; strate gy and data reduction can 
be found in lFritz et all (|2011h . 

The PACS data reduction was performed in two stages. 
The initial processing up to Level- 1 (i.e., to the level 
where the pointed photometer time lines have been cali- 
brated) is performed in HIPE v8.0 (|Qtdl2010t) using the 
standard pipeline. To remove low-frequency noise (or 
drifts) in the arrays, residual glitches and create the fina l 
maps we use SCANAMORPHOS (v!5. iRoussel l20ll . 
The final maps have a pixel scale of 2" and 3" and a 
spatial resolution of 12.5" and 13.3" full- width half max- 
imum (FWHM) at 100 and 160 fim respectively. 

The SPIRE data were processed up to Level- 1 with 
a custom pipeline script adapted from the official 
pipeline. We apply the latest flux correction fac- 
tors from the SPIRE Instrument Control Centre to 
update the maps to the lates t calibration product 
(|SPIRE Observer's Manualll201lD . For the baseline sub- 
traction we use a custom method called BriGAdE (M. 
W. L. Smith et al., in preparation) which uses an alterna- 
tive technique for correcting temperature drifts. The fi- 
nal maps were created using the naive mapper with pixel 
sizes of 6", 8"and 12" with spatial resolution of 18.2", 
24.5", 36.0" FWHM for the 250, 350 and 500 /im maps, 
respectively. All Herschel images are shown in Figure [TJ 

In addition to the Herschel data, w e make use of the 
70 um. Spitzer MIPS (|Rieke et al.ll2004| ) map published in 
IGordon et all (j2006D . This observation covers a region of 
M31 approximately 3° x 1° in size and has an exposure 
time of ~ 40 s pixel -1 . The data was proce ssed using the 
MIPS Data Analysis Tool, version 2.90 ([Gordon et al.l 
20051) and ful l detai ls of the reduction can be found in 



Gordon et all (2006). 



2.2. Gas Measurements 

To investigate the atomic hydrogen in A ndromeda we 
use th e Hi moment-zero map presented in iBraun et al.l 
(2009). The observations were made with the Wester- 
bork Synthesis Radio Telescope (WRST) covering a re- 
gion of ~ 6° x 3.5° with a resolution of 18" x 15". In 
this work we present our results using a map which has 
not been corrected for opacity effects since this correc- 
tion is uncertain. For the results which make use of the 
Hi map we also test how our results are affected by using 
a Hi map c orrected for self-opa city using the prescription 
outlined in IBraun et al.l (|2009| ). The uncorrected Hi map 
has a sensitivity of 4.2 x 10 18 cm~ 2 . 

For the molecular hydrogen gas con tent we use 
CO( J =l-0) observations presented in INieten et al.l 
(|2006T) made with the IRAM 30m telescope. This covers 
an area of 2° x 0.5° with a sensitivity of ~0.35Kkms _1 . 

All maps other than the Herschel images used for this 
analysis are shown in Figure [3J 

3. THE FIR-SUBMM SPECTRAL ENERGY DISTRIBUTION 

To investigate how the dust properties vary across M31 
we undertake a pixel-by-pixel dust analysis, using the 
Spitzer 70 £tm, Herschel ! PACS and SPIRE images. We 
first convolve the data to the same resolution as the 
500 (im map (our largest FWHM) by using a kernel which 
matches the point-spread function (PSF) in a particular 



band to th e 500 ami band. Th is procedure is described 
in detail in lBendo et al.l (|2012ft . The images are then re- 
binned into 36" pixels, chosen to be about the same size 
as the 500 fim beam so that each pixel is approximately 
independent from each other. For each band we sub- 
tract a background value for the whole map, estimated 
from regions around the galaxy. The uncertainties on 
the flux in each pixel are found by measuring the stan- 
dard deviation of the pixels in these background regions 
and adding this in quadrature with the calibration un- 
certainty. The flux errors in the majority of pixels are 
dominated by the calibrat ion uncertainty whic h we take 
as 7% for Spitzer 70 p (IGordon et al J 12001. 10% for 
PACS (see|Fritz_eialJ[2011J) and 7% for SPIRE (see Sec- 
tion [3J] for more details). 

3.1. S ED fitting 

For each pixel we fit the far-infrared-sub-millimctrc 
SED with a onc-tcmpcraturc modified black-body model 
described by: 

_ Kv M d B(u,T d ) 

° v - 752 w 

where M d is the dust mass with dust temperature T d , 
B(v,Td) is the Planck function and D is the distance 
to the galaxy. k„ is the dust absorption coefficient, de- 
scribed by a power law with dust emissivity index j3 such 
that k v oc . We assume a typical value for the ISM 
of K 3 50 M m = 0.192 m 2 kg" 1 ([Draind[200l . While the 
absolute value of n v is uncertain, its value will not af- 
fect trends with other parameters as it is a fixed scaling 
constant. The distance to Andromeda was taken in t his 
work to be D = 0.785 Mpc ([McConnachie et al.|[2005T) . 

We initially used a fixed value of f3 across the whole 
galaxy, but we found that with a fixed value it was im- 
possible to adequately fit the SEDs; we therefore let f3 
be a free parameter. To ensure the simplex fitting rou- 
tine did not get stuck in a local minimum, we ran the 
SED-fitter in two ways: first with all three parameters 
(Md, Trf, f3) free to vary; second by fixing the value of 
(3 while allowing M^ and T d to vary and repeating the 
process for all values of /3 in the range 0.20 to 5.90 in 0.01 
intervals, selecting the result with the lowest % 2 . Both 
methods gave consistent results but we created the final 
maps of the dust properties by choosing the result with 
the lowest \ 2 f° r each pixel. 

The SPIRE calibration has an overall systematic un- 
certainty of 5% due to the uncertainty in the prime cal- 
ibrator Neptune, and a statistical uncertainty of 2% de- 
termined from instrumental reproducibility; the SPIRE 
Observer's Manual recommends linearly adding these to 
give an overall uncertainty of 7%. To implement this in 
practice in our fitting algorithm we increased the uncor- 
rected uncertainty to give the same overall uncertainty 
when the errors are added in quadrature. This is imple- 
mented in the SED fitting by using the full covariance 
matrix in the x 2 calculation. 

We apply our own color correction to the Herschel 
maps by removing the Herschel pipeline 'K4' parame- 
ter and then convolving the SED with the filter trans- 
mission in the fitting process (for SPIRE the filters ap- 
propriate for extended sources are used). This method 
takes full account of all the wavelength-dependant ef- 
fects associated with PACS and SPIRE. In previous work 
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Figure 1. Herschel images used in the analysis of this paper. The images have dimensions of approximately 4.5° X 1.3°, with a tick spacing 
of 30' and centered on 10 h 43 m 02 s , +41°17'42". These images are all at their original resolution. 



HELGA II: Dust and Gas in Andromeda 



5 




Figure 2. Ancillary images for M31. The scale is the same as for the Herschel images presented in Figure [T] From top: Spitzer M IPS 
70 fim, the star formation rate (presented in G. P. Ford et al. in preparation), Spitzer IRAC 3.6 (im (presented in[Barmby et al. 200^), Hi 
and CO images as in Figure [T] The CO map (used as a tracer of H2) only covers an area of 2° X 0.5°. These images are all at their original 
resolution. 
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Figure 3. Distribution of \ 2 values from pixels fitted with the 
modified-blackbody model. The red line represents the 10% signif- 
icance value for one degree of freedom. 



(SmitheLaL|[2010) we found that there is a significant 
contribution from a warmer component of dust at wave- 
lengths < 70 /im and so the Spitzer flux at 70 /xm is 
treated as an upper limit in the fitting (i.e, if the flux 
value is higher than the model it does not contribute to 
\ 2 )- The omission of the warm component in the fit- 
ting process has a negligible effect on the derived dust 
mass as the cold compo nent dominates the total dust 
mass. iBendo et al.l (|2012l ) suggest that a warmer compo- 
nent could influence the dust emission up to wavelengths 
of 160 /im; to test this we repeated the SED-fitting by 
treating the Spitzer and PACS fluxes as upper limits and 
found it made a negligible difference to our results. 

To estimate the uncertainties on the results of our fits 
we use a monte-carlo technique. For each pixel we create 
a set of 1000 artificial SEDs, created by taking the orig- 
inal flux densities and adding a random value selected 
from a normal distribution with a mean of zero and a 
standard deviation equal to the uncertainty in the mea- 
sured flux (the correlations between the calibration un- 
certainties for SPIRE are also included). We estimate 
the error in the derived parameters for each pixel from 
the 1000 fits. We find that for each pixel there is an 
uncertainty of 20% in our estimate of the surface den- 
sity of the dust, of ±1.4 K in our estimate of the dust 
temperature, and ±0.31 in our estimate of /?. 

3.2. Results of the Fits 

In producing the final dust mass, temperature and (3 
images in this work, we only used pixels where the fluxes 
in all six bands (five Herschel & MIPS 70 /im) have a 
signal-to-noise greater than 5a. While this potentially 
causes us to miss the very coldest dust due to weak emis- 
sion in the shortest wavelength bands, we choose it as a 
conservative approach to ensure we have accurate esti- 
mates of temperature. In practice it is the 100 /im map, 
which has the lowest sensitivity, which limits the number 
of pixels in our selection. Despite this very conservative 
selection, there are still ~ 4000 pixels in our resultant 
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Figure 4. Examples of SED fits for pixels in different regions. 
The 70 fim point (blue) is used as an upper limit and the peak of 
distribution is shown by the dashed green line. 

maps. To see if our model is a statistically reasonable 
fit to the data, we created a histogram of the \ 2 val- 
ues for all pixels (Figure [3]) . As t he 70 fim flux is used 
as an upper limit (see Section |3.1[) and is usually higher 
than the best-fit model it does not usually contribute 
to x 2 and therefore we only have 1 degree of freedom 
(5data points - 3 pa rameters - !)■ The 10% significance level 
for \ 2 occurs at 2.71 which is shown by the red line in 
Figure [3] We find 9.8% of our fits have \ 2 above this 
level showing our model is an adequate representation of 
the data. We have also checked to see if our radial gra- 
dients in temperature, j3, and dust surface density are 
affected by lowering the criteria to include fluxes greater 
than 3cr and find no significant changes. 

Recent results from the Key Insight on Nearby Galax- 
ies: A Far-I nfrared Survey wi th Herschel (KINGFISH) 
presented in iDale et all (|2012f) have suggested that the 
one-temperature modified black-bo dy model underesti - 
mates the dust mass compared to the iPraine 
prescription. They attribute this difference to the con- 
tribution of warm dust emitting at shorter wavelengths. 
For our analysis this does not appear to be the case. 
First when we set the PACS fluxes (< 160 fim) as upper 
limits there is no significant change in our results. Sec- 
ond, if multiple temperature components were present 
this would bias our j3 values to lower values (since each 
temperature compo nent peaks at a different wavelength, 
iShettv et al1l2009bl) . but we mostly find higher (3 val- 
ues than expected (see Section |4}. Third the \ 2 analysis 
sugg ests the model i s cons istent with the data. In addi- 
tion, iMentuch etai] (|2012h . using a similar analysis, have 
found in M51 that the dust mass distribution is similar 
when using the lDraine fc Lil (|2007l ) or the single modified 
blackbody prescription. 

We applied the same fitting technique as outlined in 
Section 13.11 to the global flux densities, measured in an 
elliptical aperture with a semi-major axis of 2.0° and 
semi-minor axis of 0.73°. This produces a total dust mass 
of 10 7 ' 86±0 09 M Q with a dust temperature of 16.1 ± 1.1 K 
and (3 = 1.9 ± 0.3. The total dust mass from summing 
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all the pixels in our pixel-by-pixel analysis gives a value 
10 7A6 Mq, a factor of ~2.5 lower. This is expected as 
the fraction of 500 /zm flux in the pixels which satisfy 
our signal to noise criterion is approximately half the 
global value. Combining our global dust mass with gas 
measurements, we find a global gas-to-dust mass ratio 
of 72. The global temperature is consistent with these 
for other spiral galaxies obt ained using similar meth- 
ods (e.g.. [Davies et all 120111 ). The pixel- by-pixel anal- 
ysis shows a large range of temperatures and f3 values 
as discussed in detail in Section 21 Examples of fits for 
individual pixels are shown in Figure |4] 

3.3. Simulations 

To help understand the significance of our results and 
any potential biases or degeneracies in the parameters, 
we ran a Monte-Carlo simulation. Assuming the dust 
emits as a single-component modified blackbody, we gen- 
erated synthetic flux values for a range of temperatures 
and f3 values, with the same wavelengths as our real data. 
Noise was then added to the simulated fluxes with a value 
equal to the errors in the real fluxes (excluding the cal- 
ibration errors) for 2000 repetitions per T, j3 combina- 
tion. The calibration error was not included as it would 
systematically shift the fluxes for all pixels. We chose an 
input mass surface density of 0.5 M Q pc -2 as this roughly 
corresponds to the values found in the 10 kpc ring. 

The quantity that is most important for our work is 
the dust mass, which from Equation [T] is just a multi- 
plicative term. Figure [5] shows the mean and median 
mass returned for the 2000 repetitions as the input tem- 
perature is varied between 12 and 30 K for a (3 of 2; the 
error bars show the error on the mean. Between 15 and 
30 K, the mean and median dust mass returned matches 
within the errors the input dust surface density of 0.50 
M Q pc- 2 . 

At dust temperatures of 15 K and below, there are large 
errors on the returned mass, which is due to the fluxes at 
the PACS wavelengths not reaching the required sensitiv- 
ity to be included in the fit. For the actual data, we only 
included pixels in which there is at least a 5<r detection in 
all bands, which will avoid the wildly incorrect estimates 
of the dust mass seen in the simulation. To fully estimate 
the dust mass and temperature of very cold dust (T < 
15 K), we need flux measurements at >500/im. Never- 
theless, the lack of an excess at 500 /xm (Section 14. 3|) is 
circumstantial evidence that Andromeda does not con- 
tain very cold dust. We investigated if the results in 
Figure [5] were different for a different choice of j3 but 
found no systematic differences. 

By plotting the difference in the resultant temperature 
in each pixel compared to the input model temperature, 
we find that between 15 K and 25 K, the temperature 
uncertainty is *~0.6 K. Above 25 K the uncertainty in- 
creases, although in the simulation we did not include a 
70 /im point which would likely provide a constraint to 
our fits if the dust temperature was >25K. The simula- 
tion suggests that for input (3 between 1.5 and 2.4, the 
uncertainty in the returned value of (3 for each pixel is 
~0.1. As expected these uncertainties are lower than re- 
turned by the monte-carlo technique in Section [3. II as we 
have not included a calibration uncertainty. 

A degeneracy is known to exist between temperature 
and /3; Figure [6] shows that if there is an error in one 
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Figure 5. The range of mass surface densities returned from the 
SED fitting technique vers us the input temperature for the sim- 
ulated data in Section 13.31 with (5 = 2 modified-blackbody. The 
mean of the returned masses is shown by the blue points and the 
median is shown by the black points. The input mass surface den- 
sity of 0.50 Mq pc - 2 is shown by the red dashed line. To avoid 
overlapping data points the blue have been shifted by -0.1K and 
the black by +0.1 K. 

parameter this is anti-correlated with the error in the 
other parameter. As the distribution is centered on the 
origin, there is no systematic offset in the returned value 
of T or (3. As there is no systematic offset, the error in 
the mean values over many pixels will therefore be much 
smaller than the error for a single pixel. 

This simulation is based on the dust emission arising 
from a one-component modified blackbody. Fitting a 
one-component modified blackbody to pixels for which 
there are multiple temperature component s would pro- 
duce a bias towards smaller values of /3 (jShettv et all 
I2009bl) . although we would hope to detect this by find- 
ing high x 2 values. In Section [531 we show that different 
regions of M31 have different f3 — T relations and discuss 
why this is unlikely to be due to multiple temperature 
components. 

To summarise, while there is a f3-T degeneracy from 
the fitting algorithm this does not create any systematic 
offsets in the value returned. If the dust temperature 
falls below 15 K we are unable to constrain the SED due 
to the lack of data points beyond 500 /im. 

4. RESULTS 

4.1. Spatial Distribution of Dust Mass, Temperature 
and Emissivity Index 

By fitting SEDs to every pixel, we have created maps 
of dust surface density, temperature, (3 and gas-to-dust 
ratio which are shown in Figure [7] (for details on how the 
gas surface density is calculated see Section I4.2[) . The 
dust surface density distribution, un-surprisingly, is more 
similar to the maps of gas and star-formation rate than 
to the 3.6 fxm image (See Figure [TJ& HJ): which traces the 
stellar mass distribution. 

The (3 and temperature maps show significant radial 
variations. Figure [8] shows how the dust column den- 
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Figure 6. A density plot showing the correlated uncertainties be- 
tween and temperature. The uncertainties are taken from the 
simulated data with T = 17 K and (3 = 1.8. A clear anti-correlation 
is observed with the distribution centered on the correct values. To 
fully populate this graph we increased the simulation to 20,000 rep- 
etitions for this f), T combination. 

sity, temperature and (3 vary with radius (the physical 
radius is c alculated assumi ng an inclination of 77° and 
PA of 38°. lFritz et~aT]|201lD . In the star-forming ring at 
lOkpc, (3 has an average value of ~1.8 but increases with 
decreasing radius reaching a maximum of ~2.5 at a ra- 
dius of ^3 kpc. This is higher than found in global stud- 
ies of galaxies (iPlanck Collaboration|[2011aHDavies et al.l 
120111 : iDunne et al.ll201ll ). Howe ver similarly high v alues 
have been reported recently in iFovle et all (|2012D and 
iBracco et al.l ([201 If) for dust within galaxies. The value 
for the ring agrees well with early results from Planck 
([Planck Collab oration! 1201 laf) for dust in the galactic 
disk and the solar neighbourhood. The 10 kpc ring 
has an average dust surface density of ~0.6M e pc 
with a dust temperature of 18 K. Towards the very cen- 
ter of the galaxy the dust column density decreases to 
~O.O4Af0 pc -2 , j3 values fall to ~1.9 and the dust tem- 
perature increases to ^30 K. 

In Figure [8] we see a clear break in the radial variation 
for both temperature and (3 at a radius ^3 kpc. We fit 
a model with two straight lines and a transition radius 
(using a simplex routine) to the f3 results. The same 
method is used with the temperature values but we set 
the transition radius to the value obtained from the fit 
to the (3 values, which occurs at 3.1 kpc (shown by the 
dashed green line in Figure [5] or black ellipse in Figure 
[7]). At radii smaller than the transition radius, the tem- 
perature decreases with radius from ^30 K in the center 
to ^17K, with an associated increase in (3 from ~1.8 
to ^2.5. At radii larger than the transition radius the 
dust temperature slowly increases with radius while f3 
decreases with radius to ~1.7 at 13 kpc. The best-fit re- 
lationships between Td, (3 and R are shown in Figure [5] 
and listed below: 

(3= 0.15i?(kpc) + 1.98 R < 3.1 kpc (2) 

= -0.08i?(kpc) + 2.70 3.1 < R < 20 kpc 

(3) 



T d = -3.24i?(kpc) + 27.56 R < 3.1 kpc (4) 
= 0.12.R(kpc) + 16.40 3.1 < R < 20 kpc (5) 

We discuss the cause of the temperature and (3 variations 
in Section[SJ For Sections 15 . 1 1 fc I5~2l we consider the inner 
(R < 3.1 kpc) and outer regions separately. 

4.2. Radial Distribution of Gas and Dust 

In Figure [9] the radial variations of the atomic gas, 
the molecular gas and the gas-to-dust ratio are shown. 
We assume an X-factor of 1.9 x 10 20 cm -2 [Kkms -1 ] -1 
([Strong fc Mattoxlll996l) to convert a CO line flux to a 
H2 column density, although this value is notoriously 
unc ertain and has been found to vary with metallic- 
ity ([Strong et al.l 120041 : Irsraell 120051 ) . For our analysis 
of M31, the choice of X-factor makes very little differ- 
ence as the total gas is dominated by the atomic phase. 
Out of the 3974 pixels plotted in the gas-to-dust figure, 
only 101 have a molecular fraction of >50% and globally 
the molecular gas only constitutes ~7% of the atmoic 
gas (jNieten et al. 2006). To estimate the total gas sur- 
face density we include the contribution of the atomic or 
molecular gas in each pixel if the value is greater than 3c 
in their respective maps. Only 86 of our 5<r dust pixels 
are not covered by the CO(J=1-0) observations; these 
pixels are in the outskirts of the galaxy and we assume 
the contribution of molecular gas is negligible. We find 
a tight relation between the gas-to-dust ratio and radius 
(Spearman Rank Coefficient of 0.80) as shown in Figure 
[9l which is described by an exponential profile, shown by 
the red line. 

The gas-to-dust ratio increases exponentially from low 
values of ^20 in the center of the galaxy to ~110 at 
15 kp c typical of the Milky Wa y in the local environment 
(see iDevereux fc Youngl 119901 and references therein). 
Note that the values are not corrected for helium in 
the ISM. Metallicity gradients have been estimated from 
oxygen line r a tios [ 0111] /H/3, [OII]/[OIII] and R 23 by 
iGalarza et al.1 (jl999f ). They infer a radially decreasing 
metallicity gradi ent of —0.06 ± 0.03 d ex kpc -1 from the 
R23 parameter. iTrundle et al.l ([20021 ) calculate oxygen 
abun dance gradients b ased on a set of 11 Hn regions 
from iBlair et "all (|1982| ) and find gradients in the range 
of -0.027 to -0.013 dex kpc -1 depending on the calibra- 
tion used. If a constant fraction of t he metals in th e 
ISM is incorporated into dust grains (Ed mun dl l200l . 
one would expect the gas-to-dust gradient to be — lx 
metallicity gradient. We find a gas-to-dust gradient of 
0.0496 ± 0.0005 dex kpc -1 consist ent within the uncer - 
tainties to the gradient measured bv lGalarza et al.l (1999) 
(if the Hi map corrected for self-opacity is used the gra- 
dient slightly increases to 0.0566 ± 0.0007 dex kpc -1 ). 
This supports the claim that gas can be well traced 
by dust at a constant metallicity. To see if the uncer- 
tainty in the choice of X-factor could affect this result, 
we carried out the same procedure but limited the anal- 
ysis to pixels where the molecular hydrogen contribu- 
tion is less than 10% of the total gas mass. This pro- 
duced only a small change in the gas-to-dust gradient to 
0.0550 ± 0.0007 dex kpc -1 which is still consistent with 
the metallicity gradient. 

4.3. 500 [im Excess 
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Figure 7. The distribution of dust surface density, temperature and /3 obtained from the SED fitting technique and the distribution of the 
gas-to-dust ratio. The temperature and images include a black ellipse showing a radius of 3.1 kpc. The ticks are plotted at 30' intervals. 



Searches for a long-wavelength sub-mm excess (i.e. 
> 500 Mm) has been made in both the Milky Way 
(jParadis et al.1 120121 ) and for nearby galaxies. A sub- 
mm excess is important as it could suggest very cold 
dust is present (<15K) which would dominate the dust 
mass in galaxies. Other possible explanations of an 
excess include va riations in the dust emissivity index 
with wavelength (IWright et al.l H99ll: IReach et all 119951: 
IParadis et al.N2012D . different dust populations or con- 



tamination from a synchrotron radio source. Sub-mm 
excesses have been reported fr om observations of low- 
metallicity dwarf galaxies (see lO'Halloran et al.l 120101 : 
iGrossi et al.ll2"010uDale et al.ll2012D a nd spiral galaxies 
(see lZhu et al.ll2009l; iBendo et al.ll2006t) . Most detections 
have been made by combining FIR data with ground- 
based data at 850 /im or 1 mm data rather than from 
data only at <500/im. 
We searched for a sub-millimetre excess in M31 by 
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Figure 8. Results from the SED fits for each pixel plotted versus 
radius. The das hed red lines represent the best-fit linear model 
(see Section I5.1|l , the dashed green line represents the transition 
radius (3.1 kpc) found when fitting the f3 results. 

comparing the 500 /jm flux to our best-ht models. A 
500 excess is defined to be any observed 500 /jm flux 
that is greater than the expected flux at 500 fim from our 
modified blackbody fit to the data. Figure [TU] is a his- 
togram of the ratio of the excess at 500 (im and the noise 
on the 500 /im map. The distribution of the excess is con- 
sistent with a Gaussian function with a mean of 0.54<r 
and standard deviation of 0.31 a. The fact the distribu- 
tion of the histogram is a Gaussian suggests that what 
we are seeing is random noise with a fixed offset (not 
centered on zero). Two non-astronomical sources could 
explain a fixed offset, either an incorrect background sub- 
traction or a calibration error. The background correc- 
tion applied is quite small (0.2er), and thus an error in 
this is unlikely to be the entire explanation of the offset. 
The distribution of the excess is consistent with the 2% 
random SPIRE calibration uncertainty, so both factors 
together could explain the small offset. If the excess is 
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Figure 9. The distribution of atomic gas, molecular gas and gas- 
to-dust ratio versus radius. The two gas maps are plotted for all 
pixels > 3<r. The solid red line represents the best-fit exponential 
profile to the gas-to-dust ratio. 

generated from dust with a 10 K temperature and j3 of 
2, the dust mass is ~10 6 ' 58 Mq in our >5tr pixels, which 
corresponds to only 13% of the mass of the warmer dust. 
If we used a model containing dust at more then one 
temperature we would not get a useful constraint on the 
colder dust component without additional data at longer 
wavelengths. In particular observations at ^850 /im (e.g. 
with SCUBA2) would be useful to determine if any cold 
dust present. 

5. DISCUSSION 

5.1. Heating Mechanisms and Dust Distribution 

Recent studies by IBendo et all (|2010l 12012ft and 
IBoquien et al.l f|201lh have u sed Herschel colors to con- 
firm previous works (e.g.. lLonsdale Persson fc Heloul 
[T987tlWalterbos fc Greenaw alt 1996) that emission from 
dust in nearby galaxies at wavelengths longer than 
160 is mostly from dust heated by the general stel- 
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Figure 10. The value of the 500 fim excess, defined as the 500 /im 
- flux from our blackbody model, divided by the noise. The red 
line is a Gaussian fit to the histogram. 

lar population. These authors conclude that at wave- 
lengths shorter than 160 /im, the emission tends to be 
domin ated from warm e r dus t heated by newly formed 
stars. IMontalto et al.l (]2009l ) in a study of M31 using 
SpitzerMlR/FIR, UV and optical data, also concluded 
that the dust is mostly heated by an old stellar popula- 
tion (a few Gyr old) . To investigate the relation between 
our derived dust properties and the SFR and the general 
stellar p opulation, we have u sed the 3.6 /im map pre- 
sented in iBarmby et al.l (|2006l ) to trace the general stel- 
lar population (rather than the luminous newly formed 
stars that dominate the UV) and a map of star-formation 
rate (SFR). The SFR map is created from far-UV and 
24 /tin images which have been corrected for the contri- 
bution of the general stellar population to the emission 
at these wavelengths (for details see G. P. Ford et al. 
in preparation). These maps were all convolved to the 
same resolution and binned to the same pixel size as all 
the other maps in our analysis. 

We have plotted the fluxes from these maps against the 
results of our SED-fits (see Figure ITT]) ■ I n Section KTI we 
showed there is a clear break in the dust properties at a 
radius of 3.1 kpc. Therefore in Figure QTJ the pixels at a 
radius less than 3.1 kpc are shown in blue and those at a 
radius above 3.1 kpc in green. For each graph (and both 
sets of pixels) the Spearman rank coefficient is computed 
and an estimate made of its significance (see Tabled]). 
With such a large number of pixels, all but one of our 
correlations are formally significant. In the discussion 
below we have concentrated on the strongest correlations 
as measured by the Spearman correlation coefficient. 

There are strong correlations in Figure QTJ between f3 
and 3.6 /im emission. We believe these are most likely 
to be caused by radial variations in /? seen in Figure |8] 
and the decrease in 3.6 /im brightness with radius. We 
discuss the possible cause of the radial variation in f3 in 
Section O 

We find a strong correlation between dust surface- 
density and the SFR in the outer regions, but not with 



Table 1 

Spearman Correlation Coefficients for Properties of M31 



Property A Property B Region 



Spearman 
Coefficient 



P-Value 



Dust 
Surface 
Density 



3.6 (im flux 



SFR 



Inner 
Outer 
Inner 
Outer 



-0.73 
-0.16 
0.31 
0.74 



1.52 
1.53 
6.38 



x 10 
x 10 
x 10 
0.00 



-28 

-22 



Temperature 



3.6 fim flux 



SFR 



Inner 
Outer 
Inner 
Outer 



0.90 
-0.09 
0.04 
0.14 



6.17 
7.46 
5.97 
3.70 



x 10~ 59 
x 10~ 9 
x lO" 1 
x 10" 



-1!) 



3.6 fim flux 



SFR 



Inner 
Outer 
Inner 
Outer 



-0.52 
0.56 
0.16 

-0.24 



8.78 
8.99 
3.58 
2.85 



x 10~ 13 
x 10" 311 
x 10~ 2 
x 10" 



-52 



Notes. — The Spearman rank-order correlation coefficient and P- 
value (for the the null hypothesis that the two data sets are uncorrclatcd) 
for the scatter plots shown in Figure ! - ! II (values were calculated using the 
scipy.stats package and checked with IDL r_corrclatc routine). The inner 
and outer region contain 164 and 3810 pixels, respectively. 



the surface density of total stars traced by the 3.6 /im 
emission. This correlation is expected as stars are formed 
in clouds of gas and dust. In the inner region there is an 
anti-correlation between the dust surface density and the 
3.6 /im flux, which seem most likely to be explained by 
both quantities varying with radius: the 3.6 /im emission 
from the bulge increasing towards the center and the dust 
surface density decreasing towards the center. 

The strongest correlation is seen between the dust tem- 
perature and the 3.6 /im flux in the inner region. Figure 
[T2l shows a log- log graph of the two quantities and a lin- 
ear fit to the points. The gradient represents the power n 
where F 3 , 6fim cx T^ ust where n = 4.61 ± 0.15. While the 
correlation suggests the dust in the bulge is heated by the 
general stellar population, for a modified blackbody with 
f3 — 2 we would expect a gradient of 6. The difference is 
probably explained by the simplicity of our assumptions: 
that there is only a single stellar population in the bulge 
and that the bulge has a constant depth in the line of 
sight. If these assumptions are incorrect, the 3.6 /im sur- 
face brightness of M31 will only be an imperfect tracer 
of the intensity of the interstellar radiation field (ISRF). 

Looking at the temperature beyond 3.1 kpc, we find a 
weak, but still highly significant, correlation with SFR, 
suggesting that the ISRF has a significant contribution 
from star-forming regions. As most of the star forma- 
tion in M31 occurs in the 10 kpc ring this is to be ex- 
pected. This correlation can be seen in Figure [TBI where 
most but not all of the temperature peaks in the 10 kpc 
ring appear aligned with the peaks in the SFR map. In 
the same region there is a slight anti-correlation of tem- 
perature with the 3.6 /tm flux which could be explained 
by the radial decrease in the 3.6 /im flux whilst the dust 
temperature increases slightly with radius. The fact that 
dust temperature increases slightly with radius, while the 
number-density of stars, traced by the 3.6 /im, is falling 
with radius suggests that outside the bulge the dust is 
mainly heated by young stars. Nevertheless, the lack of 
a strong correlation between dust temperature and ei- 
ther SFR or 3.6 /im flux suggest the optical/UV light 
absorbed by a dust grain is from photons from a large 
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Figure 11. Scatter plots showing correlations between dust properties and 3.6 (im flux or SFR. The blue points are results at radii <3.1 kpc 
and the green data points are results at radii >3.1kpc. The Spearman rank-order coefficients for both sets of points are shown in the 
top-left corner of each plot. 

range of distances (e.g., photons from the bulge heating 
d ust in the disk). 

iBendo et al.l (|2012fi have studied FIR colour ratios in 
M81, M83 and NGC4203. They find the 250/350 /mi 
color ratio has the strongest positive correlation with 
1.6 /im emission. An increase in the 250/350 /im ratio 
would indicate ei ther an increase in d ust temperature or 
a decrease of {3. IBendo et al.l ([2012) conclude that the 
most likely explaination is the temperature effect, with 
the dust being heated by the general stellar population 
traced by the 1.6 /xm emission . We find a similar cor- 
relation to IBendo et al.l (|2012l ). only we trace the stel- 
lar radiation field with the 3.6 /im band. However, our 
SED-fitting results suggest that this is caused by a com- 
binat ion of chang e s in t emperature and /3. Note that 
since IBendo et al.l (f2012l) only use color ratios they are 
unable to discriminate between changes of temperature 
and /3. At radii greater than 3.1 kpc our results sug- 
gest the variation i n the 250/350 /mi i s mainly caused 
by a change in j3. IBendo et al.l (|2012t ) found that the 
70/160 /im color ratio has the greatest correlation with 
star-formation rate, which is evidence that there is dust 
at more than one temperature contributing to the 70- 
500 /tm emission. One possible explanation of our failure 
to find a correlation between the 3.6 /tm emission and 
dust temperature outside 3.1 kpc, instead of the negative 
correlation between 3.6 /im and f3 we find, might be if the 
Herschel emission at short wavelengths contains a con- 
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12. Log(3.6 /im) flux versus log(Temperature) for the inner 
The best-fit linear model is shown by the red line. 
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Figure 13. The color image shows the SFR image from G. P. Ford et al. (in preparation) which has been smoothed and re-gridded to 
match the maps presented in Figure[7] The contours are from the dust temperature map and drawn at 18.0, 19.5 and 21.0 K values. 
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tribution from a warmer dust component. This would 
mean our fits of a one-component modified blackbody 
would produce misleading results. However, as stated in 
Section [3 . 1 1 when we attempted the same SED fitting pro- 
cess but using all flux densities <160^m as upper limits, 
we see little difference in our results. 

5.2. Dust Emissivity and Temperature Relation 

The dust emissivity index (/?) is related to the physical 
properties of the dust grains, including the grain compo- 
sition, grain size, the nature of the absorption process 
and the equilibrium temperature of the dust. We would 
also expect to see a change in /3 due to environment from 
the processing of the grains via grain growth (e.g. coagu- 
lation, mantle accretion) or destruction through surface 
sputtering by ions/atoms or shattering by shocks. In 
M31, we detect an apparent inverse correlation between 
T(j and (3 for the inner and outer regions of M31, as 
shown in Figure [T3] We find the form of the relation is 
different for the two regions. 

Such an inverse relationship has been observed in 
the Milky Way with previous FIR- submm experiment s 
and surveys including ARCHEOPS ([Desert et al.ll2008l) . 
which showed f3 ranging from 4 to 1 with the dust tem- 
perature varying b etween 7 and 27 K, and PRONAOS 
(|Dupac et al.l 120031. which shows a variation of (3 from 
2.4-0.8 for dust tem peratures between 11 and 80 K. 
iVeneziani et al.l (|2010D used IR-mm data of Galactic high 
latitu de clouds and found a similar trend, and more re- 
cently IParadis et al.l (|2010f) with Herschel found a sim- 
ilar inverse relationship with (3 = 2.7—1.8 for Td = 
14 - 21 K for galactic longitude 5 9° (at longitude 30° 
/3 = 2.6- 1.9 for 18-23K). Recently iBracco et al.l f|2011h 
used Herschel- ATLAS observations to investigate (3 vari- 
ations in low-density, high-latitude galactic cirrus, mea- 
suring values of (3 ranging from 4.5-1.0 for 10 < Td < 
28 K. These Td — (3 relationships could be indicative of a 
problem with the temperature-/? degeneracy arising from 
the SED fitting, the presence of d ust with a range of t em- 
peratures along the line of sight (Sh ettv et al.ll2009bh . or 
real variations of the properties of the dust grains. 

On the assumption that the inverse correlations be- 
tween j3 and Td in Figure Q3] are not simply caused by 
the two variables being separately correlated with radius, 
we looked for other possible causes of the relationships: 

1. The fitting can lead to a spurious inverse co rrela- 
tion between /3 and T d (|Shettv et al.ll2009al Sec- 
tion 13 - 3[) . The most striking feature in Figure [14] 



is the clear separation in points between the in- 
ner 3.1kpc and the outer regions. To test whether 
these different distributions might be produced by 
the fitting artifact, we used the Monte-Carlo sim- 
ulations from Section 13.31 to simulate the effect of 
fitting a modified blackbody for various combina- 
tions of (3 and T d . The grey lines in Figure H5] 
show the best-fit relationships for different input 
Td and (3 combinations and clearly show that the 
two different relationships in the two regions can- 
not be obtained from a single population of dust 
grains. The green and blue data points represent 
the range of output Td and j3 for an input modified 
blackbody with T = 17.0 K (green) and T = 25.0 K 
(blue), with (3 = 2.0. A comparison of Figures H"4l 
and fl5l shows that in both regions of M31 there is a 
larger range of temperature and dust emissivity for 
the real data than that found in the Monte-Carlo 
simulation, indicating there are genuine variations 
of T d and (3 in both regions. Moreover, the fitting 
artifact cannot explain the observed relationships 
of j3 and temperature with radius (Figure [8]). 

Artificial inverse T d — /? relationships can also be 
produced if a one-component modified blackbody 
model is used to fit dust which contains a range of 
dust temperatures (jShettv et al.ll2009bT) . Since we 
are averaging through the disk of a galaxy along 
the line of sight (LOS), it is obviously possible 
that the dust contains a range of dust tempera- 
tures. While we cannot fully address this issue, 
our (3 values are higher than expected, which is the 
opposite of what happens from a LOS averaging 
of temperatures. We also find no statistical evi- 
dence from our fits that there is mor e than one 
component of dust. A lso IParadis et al.l (|2010t) and 
lAnderson et al.l (|2010| ) show inverse 1^-/3 relation- 
ships still exist in places where it is unlikely there 
is dust at more than one temperature. 

Variation of f3 with wavelength has been re- 
ported by some authors both from theoretical mod- 
els and labor atory experiments and from obser- 
vatio ns (e.g., iMenv et"al"1 120071 ; iCoupeaud et al.l 
120111 ) , with a transition around 500 /xm. In Section 
14.31 we show that there is no evidence for excess 
500 fxm emission, suggesting this is not an explana- 
tion of our results. 
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Figure 14. The variation of the dust temperature with emissivity index across M31. Data points are colour-coded for those within 
R < 3.1 kpc (blue) and those beyond this radius (green). Sol id lines show the best-fit relations for — in M31. The — relationship s 
in the literature are indicated by the dashed lines (including Dupac ct al. 2003; Desert ct al. 2008; Paradis ct al. 2010; Bracco ct al. 2011). 
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Figure 15. The variation of the dust temperature with emissivity index that ar ise f rom just the uncertainties in the measurements. The 
data shown use the simulations of the SED-fitting method, described in Section 13.31 The green and blue data points show the recovered 
values of and T<j for an input model with T = 17.0 K, = 2.0 (green) and T = 25.0 K, = 2.0 (blue). We have carried out the same 
simulation for input values of Tj and over the range T<j of 15-29 K in 2K intervals and in of 1.6-2.4 in 0.2 intervals. For each group 
of points we have fitted a line T^ oc n , which are the grey dashed lines. In these cases we have not shown the recovered values of T<j and 
0, mer ely the lines that are the best fit to the points. The red and black solid lines are the best fit models to the real data as shown in 
Figure IT41 When compared with Figure [141 it is clear that the uncertainties cannot account for the distribution in the real data. 
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To describe the Td — 8 relationships we use an em- 
pirical model of the f orm 8 = AT a , c o mmonly used 
in the literature Ce.g.. iDisert et all 120081 : iParadis et al.l 
[20Tot iPlanck Collaboration et al.ll2011h to fit the T d -B 
anti-correlation. The observed anti-correlation between 
T and 8 may arise due to a change in the physical proper- 
ties of grains including the grain optical constants chang- 
ing with temperature for amorphous grain s or changes 
in th e dust emissivity with wavelength (see lYsard et al.l 
120121 and references therein) . Other possibilities include 
grain growth or quantum mechanical effects (though 
these latter grain properties only arise at lower dust 
temperatures than observed in M31). The best-fit re- 
lationship which describes Td — f3 for R < 3.1 kpc and for 
3.1 < R < 15 kpc are (shown in Figure [T4|) : 

= f 2.30(§)- 61 R< 3.1 kpc 
\ LSB^)- 1 - 67 3.1 <R<15kpc 

where the steeper Td — 8 relationship at R> 3.1 kpc 
agrees well with the relationship found in the plane o f the 
Milky Way at longitudes of 59° (IParadis et alJl2010ft an d 
in low-density, high-latitude cirrus |Bracco ct al. 2011). 
There is some evidence that the Td — 8 relationship in 
M31 is slightly steeper, so that for the same temperature 
compared to the galactic plane, M31 has a lower 8 (but 
this is only at ~ 5 — 10% level). 

What could be a physical (or chemical) explanation 
of the different B — Td relationships in the two re- 
gions? Typical values of 8 are in the range 1.5-2.0 
for interstellar dust grains, and h ave been found in - 
global extra-gala ctic studies (e.g., ISkibba et all 120111 : 
Smith ct al. 20H lDunnelt^l201ll ) and average global 
values measured in the Milky Way (e.g., IParadis et al.l 
12015 iBracco et all [201lt IPlanck Collaboration! 12011^ 7 
Low values of 8 for large grains would typically repre- 
sent freshly-formed dust grains in circumstellar disks or 
stellar winds. Alternatively, j3 ~ 1 h as been observed in 
regio ns where small grains dominate (jSeki fc Yamamotol 
HflSnh . High values of 8 (> 2) might occur due to 
grain coagulation or to the growth of icy mantles on 
the surface of the grains in denser regio ns (|Aannestadl 
[1971 ILis et aIIIT99l IStepnik et al.ll2003h . Studies have 
also suggested high values of 8 are associa ted with very 
cold dust (T < 12 K e.g., iDesert et all 12008ft possi- 
bly caused by a change in the absorption properties 
due to quantum effects at low temperatures, increas- 
ing self-absorption in amorphous grains v i a tunnelling 
(lAgladze etatl 119961: iMennella etail 119981: iMenv et all 
2007; Par adis et a 1.1120 Hi!. 

The highest value of B is seen at the 3.1 kpc bound- 
ary between the two regions. This cannot be caused 
by changes in the quantum mechanical absorption, since 
this is only thought to be important for cold dust at 
temperatures < 12 K. The high 8 values could be due 
to efficient grain coagulation or mantle growth in dense 
molecular clouds, although this too seems unlikely as lit- 
tle CO(J = 1 — 0) is observed in this region. While 
there is no obvious explanation for the high 8 values at 
this radius, there are many indications that this 3 kpc 
"boundary" is an interesting regime, we discuss this fur- 
ther in Section l5~3l 

In the inner region, we suggest that the decrease in 8 



with corresponding increase in Td might be caused by 
the increased intensity of the ISRF. Towards the center 
of M31, we would expect increased sputtering or subli- 
mation of mantles from the increased ISRF, shown by the 
increased temperature of the dust a nd the increased X - 
ray emission observed in the center (jShirev et alJ I2001). 
The lack of gas in the central regions (Figure IH]) also sug- 
gests that dust is less likely to be shielded and thus more 
efficiently sputtered and leading to smaller grain sizes. 

As we mentioned above, a problem with this analysis 
is that it is difficult to determine which are the causal re- 
lationships. For example, we have argued that the radial 
variation in Td is due to the radial variation in the ISRF. 
The radial variation in 8 might then be due to a physical 
relationship between Td and 8 or it might be the case 
that there is no causal relationship between these param- 
eters but the radial variation in 8 is caused by a differ- 
ent effect. For example, an interaction between M32 and 
M31 might have caused a wave of star formation which 
has moved out through the galaxy, which might have led 
(by a number of processes) to the radial variations in 8. 
Therefore, we can rule out some hypotheses but we can- 
not conclusively determine which is the true explanation 
using this data set. 

5.3. Why a Transition at 3.1 kpc? 

Interpreting the transition in dust properties seen at 
3.1 kpc (Figure[8|) is difficult. One possi ble clue comes 
from previous gas kinematics studies. iChemin et al.1 
(2009) found that the Hi rotation curve inside a 4 kpc 
radius is warp e d wit h respect to the rest of disk. 
iStark fc Binnevl (|1994f ) suggest the inner Hi data is con- 
sistent wit h a bar extend ed to 3.2 kpc, while a newer 
analysis by iBermanl (|2001h explains the Hi distribution 
as the result of a triaxial rotating bulge. iBlock et al.l 
(2006), using Spitzer IRAC observations, identified a new 
inner dust ring with dimensions of 1.5 x lkpc. By using 
the stellar and gas distributions and from the presence 
of the 10 kpc ring, they conclude that an almost head- 
on collision has occurred between M31 and M32 around 
210 million years ago. This collision could explain the 
perturbation of the gas observed in the central 4 kpc. 
These other observations all show that the inner 3 kpc of 
M31 is an intriguing region, although it it not clear what 
are the causes of the difference in the dust properties. 
The perturbation of the gas may have lead to the pro- 
cessing of dust grains, or potentially material from M32 
could have been deposited after the interaction. The to- 
tal dust mass for the pixels in our selection within 3.1 kpc 
is 1O 4 2 M which is a plausible amount to be deposited 
as recently dust masses o f ~10 5 MfD have be en reported 
in Virgo dwarf ellipticals (jGrossi 

Another possibility is the dust properties are affected 
by t he differences be t ween conditions in the bulge and 
disk. iCourteau et al.l (|2011l) decomposed the luminosity 
profile of Spitzer IRAC data into a bulge, disk and halo. 
From their Figure 16 we can see that our transition ra- 
dius of ^3 kpc is approximately where the bulge emission 
begins to become a significant fraction of the optical disk 
emission. Whether the transition in the dust parameters 
is due to the changing contribution to the ISRF from the 
general stellar population and star formation or if there 
is another influence in the bulge is unknown. 
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5.4. Dark Gas and X factor 

The detection of 'dark gas' in the Milky Way was a sur- 
prising early result from Planck ([Planck Collaboration! 
I2011al) , obtained by combining IRAS 100 /Ltm data and 
the six-band Planck data from 350/xm-3mm. The 
Planck team compared the dust optical-depth with the 
total column density of hydrogen (7V^ ot ), where the op- 
tical depth at each wavelength is given by 

Tv = h (7) 

B(v, T dust ) 

where I v is the flux density in that band and B(y, Tdust) 
is the blackbody function. They assumed that at low 
Nft 0t the atomic hydrogen dominates over the molecular 
component while at high column density the molecular 
hydrogen dominates the emission. For these two regimes 
they found a constant gas-to-dust ratio, but at inter- 
mediate column densities they found an excess of dust 
compared to the gas. This excess is attributed to gas 
traced by dust but not by the usual Hi and CO lines, 
and is found to be the equivalent 28% of the atomic gas 
or 118% of the molecular gas. The excess dust emission 
was typically found around molecular clouds, suggesting 
that the most likely cause is the presence of molecular 
gas not traced by the CO line. 

We attempted t he sa me analysis as the 
iPlanck Collaboration! (|2011af ) for M31 using our 
SED fitting results from Section 2J Instead of using 
Equation we compare the column density of gas 
estimated from the Hi and CO with the column density 
of dust (for convenience we call this T,d us t)- We use 
this parameter as it is calculated with data from all 
wavelengths, whereas if we used Equation!?] small errors 
in temperature would cause large uncertainties in t v for 
wavelengths close to the peak of emission. 

The Planck team found no radial variation in the 
gas-to-dust ratio in the M ilky Way (see Figure 10 
IPlanck Collaborationl [201 lbf ). In Andromeda we show 
that the gas-to-dust ratio does vary radially (Section l4~2l 
and Figure [5]), as expected from the metallicity gradient. 
To determine if there is an excess at intermediate column 
density in dust compared to that expected from the gas 
we have to correct for the radial change in gas-to-dust 
ratio. To remove this dependence, we adjust dust col- 
umn density by using the exponential fit (shown by the 
red line in Figure so the gas-to-dust ratio at all radii 
has the same value as the center of M31. To avoid bias- 
ing this correction by assuming a value for the X-factor 
which is highly uncertain, we estimate this relationship 
from pixels where the atomic hydrogen column density 
is > 95% of Nl ot (1569 pixels out of the 3600). 

In Figure [To] we show that the relationship between 
corrected dust column density (T>dust) an d gas column 
density is well represented by assuming the two quan- 
tities are directly proportional with no excess in dust 
column-density that could be attributed to 'dark gas'. 
Although only a small proportion of the gas is molecu- 
lar, we still need to use a value for the X-factor. We can 
estimate this quantity from the data itself by finding the 
values of the X-factor and the constant of proportion- 
ality between the gas column density and the corrected 
^dust that gives the minimum \ 2 value (the fitted line 
is shown in Figure [16]) We find a best value for the 



X-factor of (1.9 ± 0.4) x 10 20 cuT 2 [K kms" 1 ]" 1 (or ex- 
pressed as aco = 4.1 ± 0.9M© pc -2 [K kms -1 ] -1 ), where 
the random error is estimated using a monte-carlo tech- 
nique (similar to one used in Section 13.11) . For the dust 
column densities in each pixel we use the uncertainties 
provided by the SED fitter as explained in Sec tion 13.11 
which is on average ~22%. INieten et al.l (|2006[) quote a 
calibration error of 15% for the CO observations (which 
directly results in at least a 15% uncertainty in the X- 
factor) which we combine with the noise in each pixel of 
our processed moment-zero CO map. For the Hi obser- 
vations we use an uncertainty map provided by Robert 
Braun (private communication) which has an average un- 
certainty of 12% on the raw 10" map. We also include 
a 5% systematic uncertainty (e.g., calibration uncertain- 
ties). If the opacity corrected Hi map is used our best fit 
X -factor is (2 ± 0.4 ) x 10 20 cm" 2 [K kms" 1 ]- 1 . 

ILerov etaLI ([2011) find values of the X-factor between 
0.97 to 4.6 x 10 20 cm _2 [K kms -1 ] -1 when analysing 
the southern, northern and inner regions for M31 with 
Spitzer data. Our average value of the X-factor falls 
within their range of values. A full analysis of the spatial 
variations of the X-factor with our data will be under- 
taken in a future paper. 

There are two main problems with this method which 
could lead us to miss 'dark gas'. First, unlike the Milky 
Way we have to average through the whole disk of M31 
and, second, M31 has a significantly lower molecular gas 
fraction. The latter prevents us from fitting the model to 
pixels with just very high and low values of Nfj ot as the 
pixels with highest molecular gas fraction are not clus- 
tered to high Njj values (this is illustrated in Figure ITT]). 
This suggests Andromeda may not be the best galaxy for 
this analysis as the molecular contribution to the overall 
column density is quite low. 

However, we can also try an alternative method of 
looking for 'dark gas' because we have one important 
advantage over the Planck team: we can see M31 from 
the outside. We can therefore make a map of the ratio of 
radially corrected dust column density (Edust) to gas col- 
umn density to look for regions of enhancements in this 
ratio (Figure [T8| . The image clearly shows spatial varia- 
tions which could either suggest regions of 'dark gas' or 
local variations in the metallicity or emissivity of dust. 
To distinguish between these scenarios an independent 
measurement of the 'dark gas' is required. On the out- 
skirts of molecular clouds, CO could be photodissociated 
and the carbon gas would therefore reside in C or C + 
([Wolfire et al.l 120101 ). Planned observations of the C[II] 
158 /im line in Andromeda with Herschel could then be a 
potential test for investigating whether 'dark gas' exists 
in M31. 

6. CONCLUSIONS 

In this paper we present the results of an analysis of 
dust and gas in Andromeda using new Herschel observa- 
tions from the HELGA survey. We have ^4000 indepen- 
dent pixels with observations in the range of 70-500 /im. 
We find the following results: 

1. We find that a variable dust emissivity index, /?, 
is required to adequately fit all the pixels in An- 
dromeda. When a variable j3 is used, the modi- 
fied blackbody model with a single temperature is 
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Figure 16. Radially corrected E^ust versus total column density of gas. The plot is shown using our best-fit value of the X-factor of 
(1.9 ± 0.4) X 10 20 cm - 2 [K kms -1 ] -1 . The red line represents the best fit model to the data assuming that Tidust °c Nj^ ot . The plot shows 
that, unlike the Planck data for the Milky Way (Planck Collaboration 2011a), at intermediate gas column densities we do not find an 
excess in dust column density over the best-fit model which would indicate the presence of gas not traced by the Hi and CO ('dark gas'). 
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Figure 17. Radially corrected E^ust versus column density of gas for pixels where the molecular fraction is greater than 20%. The data 
points arc colour-coded with the fraction of molecular gas compared to total gas. The Figure shows th at t he high column densities are not 
dominated by regions of molecular gas traced by the CO. The red line is the fitted model from Figure \TE\ 



found to be a statistically reasonable fit to the data 
in the range 100-500 /im. There is no significant 
evidence of an excess of dust emission at 500 fim 
above our model. 

There are two distinct regions with different dust 
properties, with a transition at R — 3.1 kpc. In the 
center of Andromeda, the temperature peaks with 
a value of ~30 K and a (3 of ~ 1.9. The temperature 
then declines radially to a value of ^17 K at 3.1 kpc 
with a corresponding increase in (3 to ~2.5. At radii 



larger than 3.1 kpc (5 declines but only with a small 
associated increase in temperature. 

3. The drop in (3 towards the center of the galaxy may 
be caused by increased sputtering or sublimation of 
mantles from an increased ISRF. The origin of the 
high j3 values at 3.1 kpc from the center is less clear 
but may be indicative of either grain coagulation 
or an increase in the growth of icy mantles. 

4. The dust surface density for our pixels in which 



18 



Smith, Eales, Gomez et al. 




Figure 18. Map of ^dust/Njj 01 ratio in M31. Higher values represent areas where there is less gas than predicted from dust measurements. 
The tick spacing represents 30'. 



flux is detected at > 5a in all Herschel bands range 
is between ~0.1-2.0M© pc~ 2 . We find the gas-to- 
dust ratio increases exponentially with radius. The 
gradient matches that predicted from the metallic- 
ity gradient assuming a constant fraction of metals 
are included into dust grains. The dust surface 
density is correlated with the star formation rate 
rather than with the stellar surface density. 

5. In the inner 3.1 kpc the dust temperature is corre- 
lated with the 3.6 /an flux. This suggests the heat- 
ing of the dust in the bulge is dominated by the 
general stellar population. Beyond 3.1 kpc there is 
a weak correlation between dust temperature and 
the star formation rate. 

6. We find no evidence for 'dark gas', using a sim- 
ilar technique as the Planck team. However, we 
find this technique may not be as effective for M31 
due to poor angular resolution and line of sight 
effects. We have used an alternative technique 
by constructing a gas-to-dust map after correcting 
for the radial gradient. We do find regions with 
enhancements (i.e., higher values of Sdust/A^ *), 
which may show places where 'dark gas' exists, or 
may be due to local variations in the gas-to-dust 
ratio. A detection of a potential component of CO- 
free molecular gas will be possible with future ob- 
servations to measure the C[II] line planned with 
Herschel. 

7. By minimising the scatter between our radially 
corrected dust column-density and the column- 
density of gas inferred from the Hi and CO line 
we find a value for the X-factor of (1.9 ± 0.4) x 
10 20 cm _2 [K kms -1 ] -1 (or expressed as aco — 
4.1 ±0.9M Q pc" 2 [Kkms- 1 ]- 1 ). 



Our results of Andromeda represent the largest re- 
solved analysis of dust and gas in a single galaxy with 
Herschel. The results of this analysis on M31 is strik- 
ingly different from those obtained by the Planck team 
in the Milky Way, since we find no clear evidence for 
'dark gas', a radial gradient in the gas-to-dust ratio and 
evidence for radial variation in the dust emissivity index 
(/3). In future work it will be important to understand 



these differences between the two big spirals in the local 
group. 
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